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Abstract 

We consider the k=\ Friedman-Robertson- Walker (FRW) model within loop quantum cosmol- 
ogy, paying special attention to the existence of an ambiguity in the quantization process. In 
spatially non-flat anisotropic models such as Bianchi II and IX, the standard method of defining 
the curvature through closed holonomies is not admissible. Instead, one has to implement the 
quantum constraints by approximating the connection via open holonomies. In the case of flat 
A;=0 FRW and Bianchi I models, these two quantization methods coincide, but in the case of the 
closed k=l FRW model they might yield different quantum theories. In this manuscript we explore 
these two quantizations and the different effective descriptions they provide of the bouncing cyclic 
universe. In particular, as we show in detail, the most dramatic difference is that in the theory 
defined by the new quantization method, there is not one, but two different bounces through which 
the cyclic universe alternates. We show that for a 'large' universe, these two bounces are very 
similar and, therefore, practically indistinguishable, approaching the dynamics of the holonomy 
based quantum theory. 
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I. INTRODUCTION 



Loop quantum cosmology (LQC) has become in the past years an interesting candidate 
for a quantum description of the early universe via homogeneous cosmological models [IH3] ■ 
Based on the same quantization methods of loop quantum gravity [1] , it has also become a 
testing ground for different conceptual and technical issues that arise in the full theory. It 
is perhaps not surprising that the model was first fully understood is the spatially flat k=0 
FRW cosmological model coupled to the simplest kind of matter, namely a mass-less scalar 
field that serves as an internal time parameter [5HT0]. It was shown numerically that the 
big bang singularity is replaced by a quantum bounce [7] , that connects an early contraction 
phase of the universe with the current state of expansion. By means of an exact solvable 
model, this bounce was then understood to be generic and present for all states of the theory, 
and the energy density was shown to be absolutely bounded by a critical density pcrit of the 
order of the Planck density [S] . It was then shown that semiclassical states after the bounce 
have to come from states that were also semiclassical well before the bounce [9l [TH |T2] . This 
results have also benefited from uniqueness results that warranties the physical consistency 
of the theory |T3]. The same quantization methods were applied to other isotropic models 
with and without a cosmological constant. Thus, a closed k=l was extensively studied in 
[H] and [15], while the open k=-l was considered in [16]. A detailed study of singularity 
resolution for these models was recently completed in [T8] , extending previous results for the 
flat case [T7|. For the fiat model, a cosmological constant was included in and a massive 
scalar field in [20] , where singularity resolution was also shown to emerge as a feature of the 
theory. 

An extension of this consistent quantization method was successfully implemented for the 
simples anisotropic cosmology, namely a Bianchi I spacetime in [2lj. It was soon realized 
that, for anisotropic models with a nontrivial spatial curvature, this quantization method 
based on considering holonomies along closed loops was no longer applicable. The operator 
associated to the field strength was no longer well defined on the kinematical Hilbert space 
of the theory used so far. The proposal put forward in [22] was to consider open holonomies 
to represent the connection, and then define the curvature out of the resulting operator. As 
it turns out, this quantization method has some resemblance to the quantization procedure 
known as 'polymerization' [23]. For the quantization of Bianchi IX cosmological models, it 
was also noted that this 'connection quantization' could be successfully implemented [21], 
and the singularity could also be resolved. 

A natural issue that one would like to investigate are the physical consequences of this 
'new' loop quantization. Do we have the same qualitative behavior as in the holonomy based 
quantization? This question has been satisfactorily (but trivially) answered in some cases 
where both quantizations are available. When the spatial curvature vanishes, as is the case 
of the k=0 FRW and Bianchi I models, both quantization methods coincide [221 12S] (once 
one appropriately fixes a free parameter). It is then quite natural to ask whether the same 
feature is present in other models where the intrinsic spatial curvature in non-trivial. Is 
there an important effect that the spatial curvature carries? In this respect, the k=l FRW 
model is unique to answer this question since, (to our knowledge) it is the only such model 
for which both loop quantizations exist. 

The purpose of this paper is to explore this issue in detail. More precisely, we shall 
develop the connection based quantum theory for a k=l FRW model and explore its more 
important features by using an effective description of the dynamics. We shall then compare 
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this description with that from the standard -holonomy based- loop quantization explored 
in [m [15], where the effective description has been show to correctly capture the dynamics 
of semiclassical states [14J. Perhaps somewhat surprisingly, what we find is that in the new 
-connection based- quantum theory, the corresponding cyclic universe undergoes a series 
of bounces and recollapses, but now there are two different kind of bounces. In the cosmic 
evolution, the universe alternates between these two bounces where both the density and 
minimum volume differ. Interestingly, for universes that grow to become 'large' before the 
expansion stops, the two bounces become more similar to each other, so that for a large 
universe like ours, they become almost indistinguishable. 

The structure of this manuscript is the following: In Sec. [TTjwe recall the classical k=l 
model, introducing some new notation. In Sec. Ill we recall the effective description of 
the holonomy based quantization and explore some of its consequences. Section IV is the 
main section of the paper. In the first part, we develop the loop quantization of the model, 
and in the second part we consider its effective description. We analyze then some of its 
consequences. We end in Sec. |V] with a discussion. In the Appendix we summarize our 
conventions and the computation of closed holonomies. 



II. PRELIMINARIES: THE k=l COSMOLOGY 

The spacetimes under consideration are of the form M = S x M, where E is a topological 
three-sphere S'^. It is standard to endow S with a fiducial basis of one- forms lol^ and vectors 
'fef. The fiducial metric on S is then °qab '■= l^aljlkij, with kij the Killing-Cartan metric 
on su(2). Here, the fiducial metric "qab is the metric of a three sphere of radius oq. The 
volume of E with respect to °qab will be denoted by Vq = 2%'^ a^. We also define the quantity 
io := Vq^'^. It can be written as £o =: aao, where the quantity a := (27r^)^/^ will appear in 
many expressions.^ 

The isotropic and homogeneous connections and triads can be written in terms of the 
fiducial quantities as follows. 



^ = ^"001 ; Et = l,V^X- (1) 



Here, c is dimension-less and p has dimensions of length-squared. The metric and extrinsic 
curvature can be recovered from the pair (c, p) as follows. 



qab = ^ °qab and '-fKab = ( c - y J °qab (2) 



Note that the total volume V of the hypersurface S is given by = The Poisson 

bracket for the phase space variables {c,p) is given, as in the k=0 case by, 

{c,p} = ^^, (3) 



^ Note that these conventions foUow those of [TJ (compare to [IH])- In spite of this, several of our equations 
will be different. 
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with 7 the Barbero-Immirzi parameter. From here, one can calculate the curvature of 
the connection on S as, 



ab 



c k o i o j 

^2 



The only relevant constraint is the Hamiltonian constraint that has the form. 



H 



grav 



e^h e-' E1E] F^, - 2(1 + j')e-' E^E] 



(4) 



(5) 



where e = A/|deti?|, and K^^ is the extrinsic curvature. By means of the relation = 
+ 7-?^*, with the spin-connection compatible with the triad, we can re-express the 
second term of the Hamiltonian constraint as, 



EfE'^KLK^ 



(6) 



Here fi^j, is the curvature of the spin-connection FJ^. The advantage of this substitution is 
that for this model, this expression has a simple form. 



^^ab 



C ^ °, 1* °/ I-' 

~ ~2 ^b 
"0 



With this, the gravitational constraint can be reduced to. 



8716-72 



\p\ [(c-a)^ + 7V] 



(7) 



(8) 



It is convenient to introduce new variables [S]: /3 := c/|p|^/^ and V = p^^"^. The quantity V 
is just the volume of S and (3 is its canonically conjugate. 



{/3, V} = 47rG7 



(9) 



We can then compute the evolution equations of V and /3 in order to find interesting geo- 
metrical scalars. Then, 

V = {V, -Hgrav} = ^{^y- ^^'^') (10) 

from which we can find the standard Friedman equation using the constraint equation T-L = 
^ and "Hmatt = Vp, 



V_ 
W 



8nG 



Y2/3 ■ 



(11) 



We can now compute f3 = {f3,T-L}, 



27 



where we have used the standard definition of pressure as P 
the time evolution of the expansion parameter 6 = 3H as. 



+ AnG-fP (12) 
^^'""^ We can readily find 



dV 



9 = 47iG{p - 3P) - 



3a^ 

Y2/3 



(13) 



From Eq. (11 ) we can see that the condition for a turnaround point, namely when if = is 
that the density satisfies ptum := ym- This is the point where the Hubble parameter 
vanishes. From (12) we see that, if P > — p/3 then ^ < at the turnaround point, which 
means that there is a transition from an expanding phase (where 6' > 0) to a contracting 
phase (where 6* < 0), so it corresponds to a point of re-collapse. 



III. LOOP QUANTIZATION I: THE HOLONOMY WAY 

This section has two parts. In the first one, we recall the effective equations for the 
quantization of the k=l model as developed in Ref.fH], and explore some of its consequences 
for arbitrary matter content. In the second part we restrict our attention to the case of a 
mass-less scalar field. 



A. Effective equations for holonomy-based quantization 

The basic strategy of loop quantization is that the effects of quantum geometry are 
manifested by means of holonomies around closed loops that carry the information about 
the field strength of the connection. As is shown in detail in the Appendix, the curvature 
takes then the form, 

Apfe _ sin^/i(c-a) -sin^(/ia) ^ , ■ 

where /i = ^/}?J\p\. In terms of the new variables /3 = c|p|~"^/^ and V = 1^1^''^, it can be 
written as, 

1/2/3 

X' = [sin2(A/3 -D)- sin' D] e,/ '^ul (15) 

where we have defined D := \a/V^^^. With this form of the curvature as defined by closed 
holonomies, and neglecting the so called inverse triad corrections, one can arrive at the form 
of the effective Hamiltonian, 

-Hes = - Q^(^^2X2 ^ [sin2(A/3 - D) ~ sin^ D + {1 + ^')D'] + pV (16) 
We can now compute the equations of motion from the effective Hamiltonian as, 

V = {V, nes} = {V, P}^ = y^yM>^f^ - D) cos(A/3 - D) . 
From here, we can find the expansion as, 

9 = ^ = ^ sin(A/3 - D) cos(A/3 - D) = -|- sin 2(A/3 - D) . (17) 
V A7 2A7 

From the above equation we can see that the absolute value of expansion has an absolute 
upper limit equal to \9\ < 3/2A7. We can now compute the modified, effective Friedman 
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equation, by computing — 



9 ' 



H'= I - - J p + sin^ D - (1 + 7^)D^ 1 - ^ ^ ' p - sin^ D + (1 + 7')/^ 



SttG / p- pi 
= -I^{p-Pi) 1 

(18) 

where pi = pcrit 

[(1 + 72)L)2 _ sin2 and Pcrit = S/{8nG'y^X'^) is the critical density of the 



k = FRW modeL We can immediately note from Eq. (18) that there are two points 
where the Hubble parameter H vanishes and the Universe has a turnaround. The first one 
corresponds to the point p = pi. Note that pi, in the limit A — )■ 0, tends to pi 1— 



which is the classical value for re-collapse as given by Eq. (11). Thus, in the limit of large 
volumes one expects pi to represent the density at re-collapse. The second value for density 
where the Hubble parameter vanishes is given by p = pcrit + Pi- Note that these densities, 
where there is a turnaround, is not an universal constant for all trajectories as was the case 
for the k=0 model (for the bounce at p = Pcrit)- Instead, the quantity pi is a function of 
volume and depends on each individual trajectory. The second density for turnaround is 
bounded below by Pcrit-^ There is an alternate way of analyzing the two turnaround points. 



From the expression of the expansion (17) we can see that the Hubble parameter vanishes 
when 

sin 2(A/3 - D) = sin(A/3 - D) cos(A/3 - D) = (19) 

There are two possibilities for this. 

i) When XI3-D = vr , 

for n integer, which corresponds to p = pcrit + Pi- The other possibility is, 

ii) — D = 17171 

where m is an integer number. This corresponds to p = pi. 

In fact, these considerations suggest that we could define a new variable (3 := (5 — D/X = 
(c — a)/y/p, that would also be 'conjugate' to V {{f3,V} = AtiG^). In terms of (3 many 
expressions would simplify, and it would reduce to (3 in the k=0 case. 

In order to determine which of the turnaround points corresponds to a bounce and which 
one to a re-collapse, we need to consider the rest of the effective equations of motion, 

13 = AnG-fP 
1 



-[3sm\X/3 - D) -Ssin^ D + Dsm2{X/3 - D) + Dsm2D + {1 + ^'^)D^] (20) 
27A 

= -47rG7 [p - P2 + P] 
where 

P2 = [2(1 + ^^)D - sin 2(A/3 - D) - sin 2D] (21) 



^ Also note that since pi depends explicitly on the volume, the values it takes at the bounce and classical 
turnaround point are different, so it could happen that p = pi is actually larger than in the other root, 
and it corresponds to the bounce while p — pcrit + Pi corresponds to a re-collapse [26] . 
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The Ricci scalar is given by, 



SttCp 1 + 2 



. P- Pi 

Pcrit 



327rG'pi 1 



P-Pi 

Pcrit 



2A'rG{P - Pa) 1 - 2 



. P-Pi 

Pcrit 



^2/3 

(22) 



The time derivative of the expansion is given by, 



e = cos2(A/3-D) ( + 



with 
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P3 = P2 + 



7A J 

Pcrit D , 



3 ■ eD\ 
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1 -2 



p- pi 

Pcrit 



-12nG {p - PS + P) 



1 -2 
(23) 



P-Pi 

Pcrit 



sm2{X/3-D) = [2(l + 72)D-sin2Z}] 



3 ' ' ^ 3 

Finally, the contracted Ricci curvature appearing in Raychaudhuri equation is given by. 



RabCe = -0-^0^ = ^^Gp 



Pcrit 



+8nGpi 1 - 



P-Pi 

Pcrit 



-12nGiP-ps) 1-2 



. P - Pi 

Pcrit 



It is straightforward to show that the continuity equation p + 3H{p + P) = is also satisfied 
in this case 1181. 



Let us now determine the nature of the turnaround points. From Eq. (23) we can see 
that in case i) above, where 9 = and p = pcrit + Pi, we have then. 



9 = --/3 
7 



(24) 



Therefore, the nature of the turnaround is determined by the sign of /3. If /3 < then 
^ > and the point corresponds to a bounce. However, if /3 > then ^ < and the point 
corresponds to a re-collapse. 



For case ii), again from Eq. (23), and using 9 = and p = pi we can see that, 

9 



7 



(25) 



Therefore, if /3 < then 9 <0 and the point corresponds to a re-collapse. In the other case, 
when /3 > then 9 <0 and the point corresponds to a bounce. From this discussion, we can 
see that the nature of the turnaround points can change if, during the dynamical evolution, 
P changes sign. This phenomena has indeed been observed in certain cases |26j . 



B. Concrete example: A massless scalar 

Up until now, we have considered arbitrary matter sources. Let us now restrict our 
attention to the simplest case of a massless scalar field 0, where the density is given by 
p = 0^/2 [H]. In this case, /3 < and does not change during the dynamical evolution. 
This means that the case i) above corresponds to the bounce and case ii) to the re-collapse. 



7 



In order to find the minimum and maximum volume we can put the maximum or minimum 
density in one side of the expression of density to have, 



2^2 

max 



Pcrit 



;i+7^ 



AV2 

vmax 



sm 



ymax 



and 



mm 



Pcrit 



1 + (1 + 7 



2/3 



sm 



1/3 



(26) 



(27) 



From numerical simulations performed in Ref. ^1] and analytical considerations for the k=0 

model [12], we know that the constant of the motion determines how semiclassical the 

state is. To be precise, as one increases the value of p^, in natural Planck units, it becomes 

easier to construct semiclassical states peaked on that value of the field momenta. It is then 

natural to expect that measures in a way, how large the Universe can grow before the 

re-collapse phase starts. That is certainly true for the classical equations of motion. Since 

we expect that the classical equations are a good approximation to the effective equations 

of motion in the low density regime, the volume at with the expansion stops should coincide 

when this transition happens at low densities in Planck units. Therefore, let us assume that 
1/3 

Vmix ^ crA, which means. 



~ 2 Pcrit 



(1 + f) 



AV 

y2/3 



sm 



frnax 



2V= « - ^ 

^ ''max Pent ofi 
t'max 



(28) 



from which we can see that the maximum value of volume approaches the classical value 



3/2 



(29) 



from above. Let us now estimate the value of the bounce in the same regime, where the 
value of Pcf, is large. 



2^minP 



crit 



2^2 



l + (l + 7' 



AV 

mm 



A. 



a 



— sm 



mm 



^Knin Pcrit 



(30) 



Therefore, the volume at the bounce also approaches the k=0 value 



(31) 



crit 



from above. 

In Fig. [1] we have plotted the time evolution of three universes for three different values 
of volume Vb at the bounce. From our previous expressions we see that the higher the value 
of the volume at the bounce, the higher the field momentum and the more semiclassical 
the trajectory. Note that this can be seen from the fact that the universe grows to larger 
values as one increases p^, and the density at the bounce decreases and tends to the value 

Pcrit • 

To summarize this section, we have seen that the effective dynamics of the holonomy 
based quantization, as defined in yields a cyclic universe with a bounce at a matter 
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V=500Vp| 
V=1000Vp| 
V=4000Vp| 




FIG. 1: For three values of the volume at the bounce T4), we plot the time evolution of the volume V 
(left) and the density p (right). These correspond to the values Vt, = 500^pi ( — - line), Fb = lOOO^pi 
( line), and Vb = 4000^pi ( line). 

densities that are larger than in the flat k=0 case. In the 'large volume regime', the volume 
at which the expansion of the universe stops approaches the value given by general relativity. 
Through-out the evolution, a key geometrical scalar such as the expansion of cosmological 
observers remains absolutely bounded, and is saturated by all trajectories at the end of the 
superinfiation regime that follows the bounce. These results complement those of [TH] where 
it was shown that, within this quantization, singularity resolution in generic for a large class 
of matter. 



IV. LOOP QUANTIZATION II: THE CONNECTION WAY 

For Bianchi II and IX cosmological models, where the spatial geometry has non trivial 
curvature, it was realized that the standard method of loop quantization based on holonomies 
for closed loops, was not implementable in the Hilbert space of loop quantum cosmology. A 
new quantization prescription was put forward in [22\ and also employed in [24J. The basic 
idea is to define an operator for the connection, by means of open holonomies, from which 
one can define the curvature. In this section we shall employ this quantization procedure to 
the closed k=l FRW model. 

To be precise, we define the connection by an open holonomy, from which we arrive at 
the expression for the connection 

K = 5^ H (32) 

where /i is the length of the curve which we use to calculate the holonomy along it and 
here we take j2 = ^^/X^J\p\. Just as the previous section, we shall use the variables /3 and V 
instead of c and p. 

This section has three parts. In the first one we derive the loop quantization for this 
prescription, writing in detail the quantum equations that define the theory when the matter 
is given by a massless scalar field. This resulting formalism can then be directly compared 
to that of [H]. In the second part, we consider the effective Hamiltonian and equations of 
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motion derived from the quantum theory and analyze some of their general properties. In 
the last part we specialize in the massless scalar case where we can find explicit formulae 
for some of the relevant parameters of the solutions. 



A. Quantum Kinematics 

Let us start by recalling the classical Hamiltonian constraint, 

n.iass = -^-^ [V^' - 2V'/'aP + V'/'{1 + 7^)a2] + pV (33) 

where a = t/ao = {2^)^/^ and p = pl/2V^ + U^cf)) 

As is standard in loop quantum cosmology, the gravitational part of the kinematical 
Hilbert space where the constraints are to be implemented, is given by the so called polymer 
Hilbert space [23]. In that Hilbert space, we can choose a basis of eigenstates, 

v\v) = v\v) (34) 

which is related to the volume V as follows: V = [^f^^ M ^ith K = 272/(3^373). In 
this basis, exp iXP becomes a translation operator. 

e'^^/^\v) = \v + l) (35) 

then 

smXPlv) = -{\v + 2) -\v-2)) (36) 
The quantum gravitational part of the Hamiltonian constraint operator is: 



n 



3 



/3 



(37) 

When the matter is given by a massless scalar field the quantum Hamiltonian constraint is 



n 



3 



87rG72A2 



i/^l^ sin A/3f ^^^1/4 _ 2\aV^I^ sin A/^V^^/^ + \^a\\ + 7') 



/3 



(38) 

To define the operator V^~i, we first need to define |p|~i/2 by means of Thiemann's prescrip- 
tion and, since |p|~i/2 is well defined, then we can take its cube to define ^ 

|j9|-l/2^(^) = :i_:|^|l/3||^ + l|l/3_ |^_1|1/3|^^^) (39) 

and then 

V~'^{v) = ^f{v) (40) 
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where 



f{v)= (^y\v\\\v + l\'/^-\v-l\'/^f. (41) 



The action of the Hamiltonian constraint operator on a state is given by 

- fi^dl^iv, 0) = e^(^;; 0) (42) 

where the operator O is given by 

= + 4)r/^^^^v,(, + 4; 0) - ^^^\v{v + 2)|V3v,(, + 2; 0) 

r \/\v + 21 + W\v - 21 ^2 2/. 2\i n/SiT/ ,N 

+ ~ ^^^^ + 7')l^l'/']^(^; 0) 



(43) 



The final quantum theory has a structure very similar to that of [H]. The non-separable 
Hilbert space "Hkin of the gravitational degrees of freedom is decomposed into an uncountable 
number, label by a parameter e, of superselected sectors He, each of which is by itself, 
separable. The space of solutions can be given a Hilbert space structure if one restricts 
attention to positive frequency, with respect to the internal time 0. Thus physical solutions 
ip satisfy the Schroedinger like equation, 

-i5<^^ = \/§^ (44) 

A physical inner product can be defined on the space of solutions from which the physical 
Hilbert space can be constructed. An interesting avenue would be to perform a detailed 
analysis of the solutions of this theory, along the lines of [1^. We shall leave that for 
future work. Let us now consider the effective description associated to the quantum theory 
described in this part. 



B. Effective Equations 

It is straightforward to see that the effective Hamiltonian one obtains from the quantum 
theory of the previous part, when neglecting inverse scale factor effects (as was done in [H] 
and [E]), is 

^eff = -^^^J|;^^ [(sin A/3 -D)2 + 72^2] ^pv. (45) 

It is then straightforward to compute the corresponding effective equations of motion. In 
particular, by computing V = {V,'Hes}, we can find the expression for the expansion as 

^ = A cos A/5 (sin \f3 - D) . (46) 
A7 
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From which we can find the effective Friedman equation, 



H' = ^, cos' A/3 (sin A/3 - D)' = ^(p - pO(l " — ) , (47) 

A^7 <3 Pcrit 

where pi = pcritT^-D^ and p2 = PcritD[{l + 7^)-D — 2 sin A/3]. Let us now explore what is 
the difference in the behavior of the Universe as described by these equations, compared to 
the dynamics given by the holonomy-based quantization. The first obvious observation from 



Eq. (46) is that the universe undergoes a turnaround whenever the expansion vanishes. This 
can happen either when: a) sin A/3 = D, or b) when cos A/3 = 0. The first condition can also 
be written, by using (47), as p = pi = pcrit7^-D^, and in the limit D 1 -when the volume 



is large in Planck units- corresponds to the point of re-collapse. It is interesting to note 
that, in contrast to the other quantum theory, the expression for the point of re-collapse 
here coincides exactly with that of the classical theory (recall that in the previous case, we 
only recovered this value in the large volume /momentum limit). 

Just as we had in the previous case, we expect that the nature of the turnaround points 
(whether they correspond to a bounce or a re-collapse) will be determined only after we 
consider the rate of change of the expansion (the Hubble). The second condition above, 
namely condition b) can be written as p = pcrit + P2, or alternatively, as cos A/3 = 0. Now, 
for this condition "b)", there is a crucial difference with the previous case. While in the 
effective description of the holonomy based quantization all equations were invariant under 
the mapping /3 — )■ /3 + vr/A (and therefore implementing an effective periodicity of /3 with 
period vr/A), this is no longer the case here. Even when the zeros of the term cos A/3 have 
that periodicity, the term sin A/3 — Z) does not. Therefore, there are two kind of roots for the 
equation cos A/3 = 0. The first root 'b.l' occurs when /3„ = ^^^^x^, where sinA/3„ = 1. The 
other root 'b.2' is when /3m = i^^^^t^l^^ in which case sinA/3m = —1. The important thing 
here to notice is that the density (and therefore, volume) are different in these two cases, 
which implies that there are two different kind of turnarounds of type 'b ) '. 

In order to identify the nature of these turnaround point, let us use the rest of the 
equations of motion, 

/3 = 47iG-fP - [3 sin^ A/3 - AD sin A/3 + (1 + 7^)^21 ^ (43) 
27A^ 

and, from the continuity equation, we get 

/3 = -47rG'7(p - P3 + P) where pa = ^^f^ [(1 + i')D - sin A/3] (49) 
Finally, we have the change of the expansion function given as 

o r\n 

^ = -/3(cos2A/3 + DsinA/3) + — cosA/3 (50) 
7 A7 

From this last equation we can then determine the identity of the turnaround points. For 
the different defined above we have, 

Case a): It is defined by sin A/3 = D, or alternatively by p = pi = pcrit7^-D^- In this case, 

e = - /3(cos2 A/3 - sin^ A/3 + D sin A/3) = - /3 cos^ A/3 (51) 
7 7 
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Thus, just as it happened in the holonomy-based quantization, when $ < this point 
corresponds to a re-collapse, while in the case that /3 > 0, this is a bounce. 
Case b): It is defined by cos A/3 = 0, or equivalently by p = pcrit[l + -D((l + 7^)-D — 2sinA/3)]. 
In this case we have two subcases, corresponding to the two roots of the equation cos A/3 = 0. 
Case b.l) This corresponds to the roots A/3„ = ^'^'^2X^'^ ' -^^^ ^ integer. In this case, sin X(3n = 1, 
so the change of the expansion in given by, 

e, = --^{l-D) (52) 

7 

Thus, we see that the nature of the turnaround depends not only on the sign of /3 but also on 
the magnitude of D. In the large volume regime, where D ^ 1, we have the same situation 
as in the holonomy-based quantization, namely that in the /3 < case, the turnaround point 
corresponds to a bounce (and in the /3 > case, to a re-collapse). The density is given then 
by, 

pI = P,,,,[{1-D)' + ^'D'] , (53) 
Let us nos consider the other root. 

Case b.2) This corresponds to the root A/3m. = ^^"^^'^^'^ for m integer. In this case, sinA/3„ = 
— 1, so the change of the expansion in given by, 

e2 = --H^ + D). (54) 

7 

We have the same situation as in the holonomy-based quantization, namely that in the /3 < 
case, the turnaround point corresponds to a bounce (and in the /3 > case, to a re-collapse). 
The density is given then by, 

pl = p,,,t[il + D)' + ^'D'] . (55) 

To summarize, instead of two turnaround points as in the holonomy-based quantization, 
this new quantization has the novel feature that there are three different turnaround points. 
In the case of large volume and for /3 < 0, they correspond to two bounces and a re-collapse. 
For extreme situations near the Planck scale and for certain matter content one might have 
different scenarios 1261. 



C. An example: A massless scalar 



Let us now consider as matter field a massless scalar field 0, for which /3 < and does 
not change sign during the dynamical evolution. Furthermore, we shall assume D < 1, in 
which case, the case a) above corresponds to the point of re-collapse, while the points b.l) 
and b.2) correspond to the two distinct bounces. The maximum value of volume is exactly 
given by, 

V„„=(^)"'pr (56) 

which is equal to the classical value for maximum volume for the FRW model with k=l. 
The equations for minimum volumes which correspond to the two different bounces are 



2V' 



Pcrit 



1/3 



.2 7'AV2 



2/3 



(57) 
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and 



Pcrit 



1- 



1/3 



.2 7'AV2 



2/3 



(58) 



In the limit of large field's momentum p^, since the volume is also large then we have 
D '^1. We can write the density at the two bounces as follows, 



Perit[(l + /^)' + f^'] and pl = p^,^,[{l-Df + ^^D'] 



from which it follows that, in the limit D they both tend to pcrit from above. Therefore 
the density at the bounce for both approaches with different quantization in this limit 
approaches pcrit the critical density for the /c=0 FRW model. Since both bounce densities 
have the same limit, then the minimum value of the volume for both cases goes to 



2p, 



(59) 



crit 



therefore, when the field's momentum is very large, since we can ignore the negative 
powers of volume, the maximum absolute value of expansion for the second approach goes 
to 3/27 which is the same as in first approach. 





4e+04 




3.5e+04 




3e+04 




2.5e+04 






E 




Z] 




2e+04 


> 






1.5e+04 




1e+04 




5e+03 



V=500Vp| 
„ V=1000Vp| 




4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 




I hi 



V=500\/p| 
V=1000Vp| 
V=4000Vp| 



'VIA \ 



KAL 



20 



40 60 
Time 



80 



100 



FIG. 2: For three values of the volume at the bounce Vh, we plot the time evolution of the volume V 
(left) and the density p (right). These correspond to the values = 500£pi ( — - line), 14) = lOOO^pi 
( line), and Vb = 4000£pi ( line). 

In Fig. (|2| we have plotted the time evolution of the universe for different values of the 
minimum volume at the bounce. As we can see, as we increase this value, and therefore, the 
field's momentum p^, the two bounces tend to each other, both in terms of the value of the 
volume and in the maximum value of the densities. Note that the densities at the 'strongest' 
bounce are much higher, in this regime, than in the holonomy-based quantization, and that 
they decrease as one increases the value of p^. One can further compare both description 
by fixing the value of and comparing the time evolution of volume and density. We have 
plotted such comparison in Fig. ([3| for = 10^. Note that the density at the bounce in 
the holonomy-based quantization is in between the two densities for the connection-based 
quantization. The period between the point of re-collapse is not the same for both schemes 
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FIG. 3: Wc plot the time evolution of the volume V (left) and the density p (right), for the two 
quantization methods, for = 10^. 

but, as one increases p<^, they approach each other, just at the volume and density at the 
bounce converge. 

Let us summarize the results this section. First, we developed the quantum theory for 
k—1 loop quantum gravity coupled to a scalar field, employing a quantization method the 

uses open holonomies to regulate the field strength appearing in the constraint. In the 
second part we derived some of the consequences of such a quantum theory, by means of 
its effective description. We found that the most dramatic difference from the quantization 
previous explored is that the cyclic universe undergoes cycles of contraction and expansion, 
but alternating between two different quantum bounces (or alternating between two kinds 
of points of re-collapse and a bounce). Furthermore, we saw that for 'large universes', where 
the universe expands to a large volume (in Planck units), the densities (and volumes) of the 
two distinct bounces approach each other and converge to the values attained in the A;=0 
theory. 



V. DISCUSSION 

In this article we have explored a quantization ambiguity that exists for certain models 
in loop quantum cosmology. This correspond to the freedom of using closed holonomies 
around loops to define curvature or open holonomies to define connections. Since it is only 
the latter choice that is available for anisotropic models with non-trivial spatial curvature, it 
is important to understand the particular features of this quantization, and compare it to the 
original holonomy-based loop quantization. In this regard, the isotropic k=l FRW model 
is ideal since both quantizations exist and are not equivalent (while they are in the case of 
k—0 and Bianchi I). We have explored some of the differences between these two theories, by 
means of their corresponding effective descriptions. The equation of motion for both theories 
are not the same, and therefore their underlying dynamics is different. The most dramatic 
difference is that, while the universe is cyclic in the holonomy-based quantization with a 
bounce followed by a re-collapse, in the new quantization the situation is more complicated, 
with three different turnaround points. In the semiclassical limit where the universe is a 
assumed to grow large, we have seen that there are two kinds of bounces with different 
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densities that alternate with the re-collapse. The volume at which the expansion stops and 
the universe starts to contract is also different. 

Interestingly, in the limit of large universes both theories converge and the two distinct 
bounces of the connection-based theory approach that of the holonomy-based quantization. 
In this limit both descriptions approximate general relativity during the small density epochs 
of the cyclic universes, making them almost indistinguishable. It would be interesting to 
explore further the similarities and differences of the two approaches regarding singularity 
resolution, as was done in for the holonomy based description. Further numerical anal- 
ysis with various matter fields might yield significant differences that could have potential 
observable consequences. This shall be reported elsewhere [26] . 
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Appendix A: The three sphere, holonomies and curvature 

For a 3-sphere with radius equal to Oq, the line element can be written as 

ds^ = al{da'^ + d/3'2 + d7'2 + 2 cos/^da'dV) 

where < a' < n , < /3' < n /2 and < 7' < 27r. With a simple redefinition of coordinates, 
a = 2a', (3 = 2/3' and 7 = 27', it can be written as 

ds^ = ^ (da^ + d/32 + d-f'^ + 2 cos /3dad7) (Al) 

where < a < 27r, < /3 < tt and < 7 < Att. For this metric, the volume of S is 
Vo = 2(j)'^ oq. Recall that we have defined io = Vo^^, and a = io/cto = (27r^)^/^. 

Let us now compute the holonomy along the edge e with length £', parameterized by i, 
tangential to vector = (d/di)"". It is given by 

/i(^) = expiJ^A ■ de{i)) = exp(^ t^'Airjdi) . (A2) 

If we want to use some angular parameters like 6 instead of i we will have, for a general 
integral, 

di t{F) = d9 t'{F) (A3) 



with t' = 4n and a playing the role of a 'radius', since i = a9. For our problem, we can 



de 
define 
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Therefore, to calculate a component of of the curvature, we can construct a closed loop 
as follows. In coordinates (a,/3,7) 

i) Move from (0,7r/2,0) to (0, 7r/2, 2a/i) following = d/dj, 

ii) Then move from (0, 7r/2, 2(T/i) to (2(T/i, 7r/2, 2(T/i) following —^3 = d/da, 

iii) Next, move from {2afi,TT /2,2cr^) to (2(j/i, 7r/2, 0) following —"e^, and finally 

iv) Move from (2a/i, 7r/2, 0) to (0,7r/2,0) following 4- 

The open holonomy along one edge, with parameter fi is given by 

h^f"^ = exp( / t'^'AiTjde) (A4) 



^0 

where 6* = a or 7 depending on the edge, and the effective radius of the 3-sphere used to 
translate from lengths to angles is ao/2 (compatible with the fiducial metric (Al)). Thus, 
we will have for the closed loop defined above, 

then we have 

"erelF^, = lim 4^Tr(/.n3,r'=) = -l(c^ - 2ac) (A6) 

recovering thus the classical expression for curvature. If we do not take the limit /z — )■ 
but instead take the the smallest eigenvalue of the area operator, or equivalently 

fL^lpl = \^ then the curvature can be approximated, at scale A, as 

Apfc sin^ /i(c - g) -sin^ (/ior) 

F., = (A7) 

where /2 = a/ A^/|p|. 
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